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Abstract 

Metallic nanoparticles embedded in stimuli-responsive polymers can be regarded as nanore¬ 
actors since their catalytic activity can be changed within wide limits: the physicochemical 
properties of the polymer network can be tuned and switched by external parameters, e.g. tem¬ 
perature or pH, and thus allows a selective control of reactant mobility and concentration close 
to the reaction site. Based on a combination of Debye’s model of diffusion through an energy 
landscape and a two-state model for the polymer, here we develop an analytical expression for 
the observed reaction rate constant /^obs- Our formula shows an exponential dependence of this 
rate on the solvation free enthalpy change AGsoi , a quantity which describes the partitioning 
of the reactant in the network versus bulk. Thus, changes in AGsoi, and not in the diffusion 
coefficient, will be the decisive factor affecting the reaction rate in most cases. A comparison 
with recent experimental data on switchable, thermosensitive nanoreactors demonstrates the 
general validity of the concept. 
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INTRODUCTION 


Metallic and oxidic nanoparticles have been the subject of intense studies in recent decades because 
of their catalytic activity.^ For example, gold becomes an active catalyst for oxidation reactions 
when divided down to the nanoscale.l^ Titania nanoparticles are also highly active catalysts^^ 
and there is a large number of other nanometric systems with promising catalytic properties.^ Use 
in catalytic reaction requires a simple and secure handling of nanoparticles by a suitable macro- 
molecular carrier system, in particular when working in solution. Such a carrier system should 
not impede the catalytic activity of the nanoparticles but keep them firmly stabilized throughout 
the reaction and the subsequent work-up. Reactants should also be able to quickly enter and leave 
the carrier. A catalytic system thus defined acts in many ways as a nanoreactor, that is, it contains 
and shelters the catalytic reaction. Examples of these nanoreactors are given, e.g., by dendrimeric 
systems!^ or those based on spherical polyelectrolyte brushes.l^ 

In recent years, a new class of carrier systems is emerging that can be termed active nanore- 
actors.fi^^ Here the nanoparticles are embedded in a polymer gel that reacts to external stimuli. 
In solution, diffusion within the gel may be manipulated by parameters such as temperature or 
pH. The best-studied examples of such active carriers are colloidal gels made from crosslinked 
poly(N-isopropylacrylamide) (PNIPAM) that undergo a volume phase transition at 32‘^C. It has 
been demonstrated that the catalytic activity of metal nanoparticles embedded in such gels can be 
manipulated using temperature as the external stimulus.^ 

In particular, a simple model system has been prepared by enclosing a single gold nanoparticles 
in a hollow PNIPAM sphere, a so-called yolk-shell architecture.^ This active carrier in aqueous 
phase provides the ideal means to investigate the manipulation of the activity of a nanoparticle in a 
stimuli-responsive nanoreactor. It has been demonstrated that the reactivity of the nanoparticle can 
be switched depending on the hydrophilicity of the reactants.^ These findings open a new way 
to introduce selectivity into the catalysis with nanoparticles. Another model system, introduced 
by Liz-Marzan and co-workers,^ is made up by a single gold particle embedded in the middle 
of a colloidal PNIPAM sphere, a core-shell architecture. For this system, Carregal-Romero et al. 
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demonstrated that the reduction of hexacyanoferrate (III) ions by sodium borohydride is slowed 
down drastically above the temperature of the volume transition, that is, when the gel has shrunken 
and most of the water has been expelled.^ These workers explained this decrease of the rate 
constant by the decrease in the reactant’s diffusion coefficient in the gel in its dense, collapsed 
state using a two-state model for the gel. In this way Carregal-Romero et al. presented the first 
theory of the kinetics of catalysis by nanoparticles in such an active nanoreactor.^ 

Studying the reactivity of nanoparticles in the condensed phase requires a model reaction that 
allows us to obtain kinetic parameters related to catalysis with the greatest possible precision. 
Based on the pioneering work of Pal et al^ and Esumi et al.^ we)22H25] many others^^^EMsS] 
have demonstrated that the reduction of nitroarenes and especially of nitrophenol to 4-aminophenol 
by borohydride ions in aqueous phase fulfills all the requirements for such a model reaction.12^ 

Using this model reaction, Wu et al. were able to show that the reactivity of the rather hy¬ 
drophobic nitrobenzol is even increased when raising the temperature above the temperature of 
the volume transition.^ This finding cannot be rationalized anymore in terms of a changed dif¬ 
fusion coefficient for the reactants. Wu et al. called attention on the fact that the thermodynamic 
interaction of the reactants with the PNIPAM network must be considered in this case, and used 
the Debye-Smoluchowski theory of diffusion on an energy landscapej ^^l^^ l to explain the observed 
results. 

Here we present the full theory of nanoreactors that combines our previous analysis^ with a 
two-state model^i^^^llSQ] takes into account the thermodynamic transition within the gel. The 
predictions of this model are compared to recent experimental data for both yolk-shell^ and core- 
shell^ systems. The remainder of the paper is organized as follows: First, in the Theory section 
we present a step-by-step derivation of our model, carefully stating its underlying assumptions. 
We then perform in the Discussion section a parameter study to highlight the effects of the vari¬ 
ous parameters in the model, and discuss the main predicted trends providing a comparison with 
available experimental data. Finally, a brief Conclusion will wrap up all results. 
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THEORY 


Kinetics of surface reactions 

We start with the rate equation describing a chemical reaction of first order: 

= -^expCo(t) (1) 

where cq is the concentration of reactants in solution, and ^exp is the experimental rate constant. 
By using a microscopic description we will arrive at an expression for the number of reactants 
transformed by one nanoreactor per unit time at a constant bulk density cq, ending up with an 
equation of the form: 


dN 

~7~ — ^obs^O- 

at 


( 2 ) 


In order to connect ^exp and ^obs> one needs to assume that nanoreactors do not interact with 
each other and can be treated independently (i.e., we can use a cell model, as explained in more 
details elsewhere!^. Given that under relevant experimental settings low nanoreactors’ concentra¬ 
tions are used, this condition should be typically satisfied. Within this approximation: 


dN 

dt 


^obs^O 
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leading to 


^exp — ,, ^nano^obs (4) 

Vsol 

where A^nano is the total number of nanoreactors present in a volume ygoi of solution. Eq. Q 
provides a link from the rate ^exp typically used to describe experiments and the microscopic rate 
at which nanoreactors transform reactants. 

Since the reaction catalyzed by the nanoparticle takes place directly on its surface, the reaction 
rate per total surface area of catalyst (i.e., total nanoparticle surface area) and per volume of solu- 
tionl^SES^SESlis necessary to compare different systems in an unbiased manner, hence: 


j _ ^exp _ iVnano^obs _ ^obs 

~ ■ 5tot ■ ^ 

where 5o is now the surface area of a single nanoparticle inside the nanoreactor. Via Eq. (??), 
we see how this latter normalized rate constant ^exp is truly an intensive measure of the properties 
of a single-nanoreactor, and is the quantity that should be used for comparison within different 
systems. We next proceed to describe our microscopic model for the calculations of ^obs- 

Model 

A schematic representation of a yolk-shell nanoreactor is shown in Fig.[T| For such a spherically 
symmetric system, we intend to model the rate for the catalytic reaction occurring at the surface 
of the nanoparticle, which requires taking into account the diffusion process that brings a reactant 
from the bulk solvent to the surface of the nanoparticle inside the nanoreactor’s core.^ The ki¬ 
netics of this process is governed by two key functions, the diffusion coefficient of the reactants 
along the trajectory expressed by a diffusivity profile D{r) and a thermodynamic free enthalpy of 
solvation profile AGgoi (r). This local free enthalpy AGgoi (r) measures the difference in the local 
excess chemical potential and the ’reservoir’ chemical potential of the reactants far away from the 
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Figure 1: Schematic representation of a yolk-shell nanoreactor and central quantities for our theory. At the 
center core of the nanoreactor sits a metal nanoparticle (blue) of radius /^np, embedded in a spherical polymer 
shell (’cage’) of inner radius Rg and outer radius Rg + d. The reactants have to overcome a solvation free 
enthalpy barrier AGsoi (r) to reach the nanoparticle depicted by the dotted lines. In our model, AGsoi (r) is 
represented as a piecewise constant function, of value AGsoi inside the polymer shell, and zero otherwise. In 
the case shown, AGsoi is higher in the gel than in the bulk solution, representing “polymer-phobic” reactants, 
with an energetic penalty to absorb inside the shell, but the opposite situation can also occur. 


nanoreactor. Along its trajectory, the reactant encounters a different chemical environment in the 
gel, and its average interaction with it then results in local D(r) and AGsoi (f) different than in bulk. 

In order to keep the model as simple as possible, we take both D{r) and AGsoi (f) to be a step 
function centered in the polymer cage of the nanoreactor, of width equal to the shell width d (see 
Fig.[2for reference). This is a good approximation for spatially homogeneous gels that have been 
under consideration recently.^ Thus, we get: 


{ Dq for r < Ra or r>Ra + d 

( 6 ) 

Dg for Rg<r<Rg + d 

and 

{ 0 for r < or r > -|- r/ 

(7) 

AGsoi for Rg < r < Rg + d 

where Rg is the inner radius of the polymer shell and d its width, cf. Fig. 1. The choice of step 
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functions basically implies that a reactant sees only two distinct environments, the interior of the 
polymer shell or the bulk solution. Hence, in this form, AGgoi can be identified with the difference 
in the solvation free enthalpy of the reactants in the gel with respect to that in the bulk solvent. As 
such, AGsoi does not depend solely on the reactant-polymer interaction, but also on the solvent con¬ 
stituting the bulk solution. By measuring the equilibrium partitioning coefficient of the molecule 
in the gel K^q, AGsoi can be measured experimentally, or computed via atomistic simulations, with 
the following equation: 

^eq = C^lcQ = exp[-j8 AGsoi ] (8) 

where Cg and cq are the equilibrium concentrations of reactant in the polymer gel and in the bulk 
solution, respectively, and j8 = \/k^T is the thermal energy. 

As for the free enthalpy, the two diffusion coefficients Dg and Dq appearing in Eq. (??) can be 
identified as the bulk value in the polymer gel and solution, respectively. The diffusion of solutes 
through hydrogels is a complex process^^lSl and depends not only on the individual properties 
of the polymer, solvent, and solutes, such as size and concentration, but also on the particular 
interactions of the solute with the polymer network. In our case, we restrict the discussion to ho¬ 
mogeneous gels with semiflexible polymers in the dilute and semi-dilute regimes with very small 
diffusing solutes with a size comparable to a polymer monomer. In this case obstruction or hydro- 
dynamic theories for small diffusing particles should better describe our problem.S^lSI Regardless, 
the important point we want to emphasize is that all theoretical approaches to the diffusion process 
in our regimes predict a decrease of diffusion for higher packing fractions. The detailed depen¬ 
dence of Dg on the system parameters is thus not of importance to discuss the qualitative features 
of the temperature dependence of the reaction rate of nanoreactors. 

Coupling spatial diffusion and reaction: Debye-Smoluchowski approach 

Once the diffusivity profile and the underlying free enthalpy on which the diffusing reactant must 
move to reach the surface of the nanoparticle are specified, it is possible to use an approach, initially 
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developed by Smoluchowskil^^land then extended by Debye,to calculate the diffusion-controlled 
part of the reaction rate in our system. This approach was first used to describe the rate of collision 
between charged ions in solution,but the underlying physical picture is equivalent to our model. 
For the reader’s convenience we provide a step-by-step derivation of the final equations in the 
Supplementary Information, and we also suggest to consult the reviews of Calef and Deutchl^ 
and Berg and von Hippell^ on diffusion-controlled reactions. The Debye-Smoluchowski model 
describes the rate at which a particle, driven by gradients in the chemical potential, diffuse from 
a bulk solution kept at constant concentration cq towards a fixed sink of radius i?np- Close to the 
sink, a certain number per unit time of the molecules arriving are allowed to absorb (’react’) as 
expressed by a surface reaction rate constant ^r. Restated in our language, the sink is nothing 
but the nanoparticle inside the nanoreactor, and absorption into the sink means reacting in the 
proximity of the nanoparticle’s surface. The total reaction rate ky, in units of reacting particles per 
unit time, arising from this model can be written as: 


kt ^ = A:r ^ ^ or 

Tt = Tr + Td (9) 


where ^r and are the so-called reaction and diffusion rates, respectively. The times T,, Tr and Tq 
defined in Eq. (??) are thus just the reciprocal of the corresponding rates, and should be interpreted 
as the effective time to diffuse to the sink Tq, and to react once in the surface proximity in Tr. 
Under these conditions the relation between the diffusivity and free enthalpy landscapes D{r) and 
AGsoiir), respectively, and the diffusion rate is given by the Debye-Smoluchowski expression 
(see the SI) 


kj) = 47tco 


■ r exp [/3AGsoi (r) ] ' 

Anp D(r)r2 


( 10 ) 
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Solving Debye’s equation for yolk-shell nanoreactors 

In an integral of AGgoi (r) over space appears, which can be analytically expressed using our 
simplified step-wise forms for D{r) and AGsoi (f) valid for nanoreactors, see Eq. (??), obtaining: 


^-1 


td 


f°° exp[j8AG(r)] 
47tr^coD{r) 


vnp 




1 


Anr^coDo 

1 


+ 




explpAGsoi] 


R„ Anr^coDg 


iR.+d Anr^coDo 




( 11 ) 


( 12 ) 


where again we split the typical time Td into three different contributions: the effective time to 
arrive from the bulk solution to the gel. Too, that to cross the gel Tg, and that to get to the surface of 
the nanoparticle once the gel has been crossed, Tnp. 

Solving for the integrals in Eq. (??) gives the following result: 

^ _ 1 Rg ^np 

“ CqATIDo RgRnp 
_ exp(j6AGsoi) d 
^ cqAkD^ (-^g d)/?g 

Too = , ^ ^ (13) 

CQAKDoiRg + d) 

For typical nanoreactors, Rg^, /?np ( which is exact when the reactor has a core-shell in¬ 
stead of a yolk-shell structure^. Considering also that Dq > Dg, one can further simplify Eq. (??) 
as: 
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^np ~ 0 


( 14 ) 


. _ exp(/3AGsoi) 

~ coAnD^Rg 


> Too 


(15) 


Since the largest typical time, i.e. the slowest rate, dictates the final effective value for the 
reaction rate, the latter inequality implies: 


~ l/'fg = 4;rDgCo/?g exp (—jSAGsoi) 

= AncoRgB^ (16) 

where we introduced the quantity ^ = DgCxp (—jSAGsoi)? sometimes referred to as “perme¬ 
ability” (or inverse diffusive resistance),l2SESl which takes into account the compounded effect of 
the diffusion coefficient and the solvation free enthalpy. 

The dependence of the surface reaction rate on temperature can be simply modelled using an 
Arrhenius form to account for the effect of the activation energy necessary for the reaction to occur. 
In this case (defined via solution of the Debye problem, see SI, Eq. (??) is simply augmented 
by an exponential term, i.e: 

kR kRCxp > ( 17 ) 


where AEr is the activation energy of the surface reaction (see the discussion of this point in ref .12^. 

Coupling the Debye-Smoluchowski approach with a two-state-model 

The stimuli-responsive polymer gels used in active nanoreactors exhibit a volume transition be¬ 
tween a swollen and a collapsed state at its lower critical solution temperature (LCST), Tq. Since 


10 






the transition can also be induced by changing pH or other external variables, for the sake of 
generality we describe the state of the gel as being a function of an unspecified general exter¬ 
nal parameter ^, and represent it via a two-state modelMore precisely, calling AGab{^ ) = 
AGa (^) — AGb ((^) the difference of free enthalpy between the swollen (A) and collapsed (B) states, 
the Boltzmann-weighted probability to be in state A is given by: 


exp[-^AGAg((^)] 

l+cxp[-(5AGABm' 

In a two-state model, the probability of the system being in the B state is thus 1 — pa(<^)- If one 
now assumes that a molecule coming from the bulk solution towards the gel will see one of the 
two environments with a probability given by Eq. (??), we arrive at the following equation for the 
total observed reaction rate: 



kohs = ik) = PAf4 + - PA)kf. (19) 

where = k^^(AGsoi,a) + ^DH'^^sol,a:f^g,a), with a =A,B. Eq. (??) is equivalent to con¬ 

sidering a situation where a reactant diffuses in a constant environment. Each environment will 
give rise to a different rate, and the observed one is just the thermodynamic average between the 
two. This assumption, also implicitly made by Carregal-Romero et al.in their model,^ becomes 
better the larger is the time it takes a system to swap between these two states of the gel compared 
with the average time a single reactant molecule takes to diffuse through the gel (which is of order 
d^/Dg, and should not be confused by Td). 
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Figure 2: Probability to find the system in the A state, Pa{^) as a function of the control parameter 
in a thermodynamic two-states model, for two different values of the transition enthalpy AHab- 
The latter is defined through its connection with the transition free enthalpy via the standard rela¬ 
tion AGab{^) = AlHab — TASab- At the critical value for the control parameter pA moves from 
0 to 1, with increasing sharpness for increasing AHab (note that is simply 1 — pa). 


On the opposite limit instead, when the gel switches infinitely fast between the two states, a 
reactant crossing the gel would see an average environment, i.e., an average solvation free enthalpy 
and an average diffusion coefficient.In this case, we have the following equation for the predicted 
rate: 


^obs,fast — ktiiAG^oi), (Dg)), (20) 

where, as for Eq. ??, the two-state average (..) of a quantity X is again simply defined via the mean 

(X) =pA(r)XA + (l -pa{T))Xb. (21) 

There is a quantitative difference in the reaction rates calculated within these two different limits, 
as shown in Fig.|^ but both show the same qualitative trends. 

We note, however, that in those cases where the swapping rate is finite (but not zero) and close 
to the crossing rate, resonance effects might give an additional signature close to the transition 
temperature (cf., a review on stochastic resonancel^. The investigation of those effects is out of 
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Figure 3: Comparison of the reaction rate of a single nanoreactor calculated in the limit of in¬ 
finitely fast (blue) and infinitely slow (red) swelling of the polymeric gel i.e. using Eq. (??) and 
Eq. (??), respectively, for two different values of the transition enthalpy AF/trans- The two limits 
describe qualitatively similar results, whereas quantitatively their agreement depends on the tran¬ 
sition enthalpy, i.e. from the sharpness of the transition between the swollen and collapsed state as 
a function of the control parameter ^ (here, ^ is taken to be temperature). 


scope of this paper and will be addressed in future studies. 

Let us now describe the reaction rate behavior resulting from Eq. (??); close to the critical 
value for the parameter i.e., where AGab{^) = 0, the system continuously but sharply changes 
from the A to the B state (see Fig.|^. These two states of the polymer cage represent different 
physicochemical environments, hence we expect the solvation free enthalpy AGsoi for a reactant 
to also change markedly between them. Since AGsoi enters exponentially in the observed reaction 
rate, Eq. (??), a sudden change upon polymer transition of AGgoi will be reflected by a concurrent 
sudden jump in kobs- However, given that also the diffusion coefficient will jump from two distinct 
values, the overall sign of the change will depend on the permeability ^ and not simply AGgoi- 
Despite this, we notice that this latter quantity enters i)^explicitly exponentially, whereas the dif¬ 
fusion coefficient only linearly, and it is thus expected to play a larger role. 

Eq. (??), together with the definitions in Eqs. (??) and (??), allows us to clarify the connection 
between our model and the model presented by Carregal-Romero et al^ If we set AGsoi = 0, 
that is, had we not accounted for the interaction between polymer and reactant, Eq. (??) would be 
exactly equivalent to their model, which describes the reaction rate purely in terms of diffusion. 
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Eq. (??) further suggests another insightful way to interpret our results. The diffusion rate A:d 
in both models can be written as ko = AKRDc{Rg), where c{Rg) is the local concentration of the 
reactant in the gel. In the model by Carregal-Romero et al, c(Rg) = cq, i.e. the bulk concentration. 
In our model, since we consistently account for the thermodynamics of the system, we obtain 
c(Rg) = coexp[-/3AGsoi ]. 
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Figure 4: Parameteric study of the influence of various factors on the reaction rate of nanoreactors 
as a function of temperature. A) For different AAGsoi = AGsoi,s — AGsoi^ (and fixed value for 
AGso 1 ,b- B) At constant AAGsoi(< 0 ), but changing the final value of AGso 1 ,b via its enthalpy 
A/fsoi,s- C) For varying Db/Da, i.e. for different drops in the diffusion coefficient in the collapsed 
state, in the case where AAGgoi is constant (and < 0). This latter case shows how a large drop in 
diffusion coefficient can mask the effect of a decreasing free enthalpy. In this case, however, the 
apparent drop in kobs will still be smaller than what would be calculated using an ideal diffusion 
model. For all but panel C), = 0.2Z)g ,4 = Dbuik- 
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DISCUSSION 


In the following, we shortly present a parametric study to demonstrate the influence of the various 
parameters in our model. For ease of discussion, and to make contact with a known and well 
studied system, that is PNIPAM-based yolk-shell nanoreactors, we take the external parameter 
through which AG^eis tuned to be temperature. Thus, the two two states of the gel are the 
low temperature swollen state A and the high temperature collapsed state B, and the variation in 
temperature of AGsoi is described via AGsoi, Of = AHa — TASa (oc = A,B), with AHa,ASa being 
constants. 

In order to highlight interesting trends dependent on AGgoi, we take the case where the re¬ 
action rate is dominated by diffusion, i.e., kt ~ ko << ^r. Fig.|^ shows kobs vs temperature 
(normalized by the transition temperature) for three different scenarios. Panel A) shows kobs? 
normalized by its maximum value within the plotted temperature interval kmax, for changing 
= AGsoi,A “ ^GsoI,b (i-c. the difference in solvation free enthalpy between the A and B 
state). When AAGsoi > 0, both terms appearing in ^ contribute to a drop in the rate, but it should 
be noticed that this drop is much stronger than it would be observed based on the drop in the dif¬ 
fusion coefficient alone. Indeed, this fact can explain the large drop of the “effective” diffusion 
coefficient necessary to fit rate-vs-temperature curves in the model of Carregal-Romero et 
which is not compatible with predictions from the hydrodynamic theories of diffusion that should 
be valid in this experimental regime.l^^M] when AAGgoi < 0, will depend on two contrasting 
effects: a decrease in the diffusion coefficient and an increase in the local concentration of reactant 
due to better solvation in the gel. In this case, for negative enough AAGsd, a jump to higher rates 
is observed close to the transition temperature. 

Whereas AAGsoi dictates whether a jump or a drop in the rate occurs at the critical temperature, the 
absolute value of kobs depends on the absolute value of AGsoi- fo Panel B, we fix AAGsoi(< 0) by 
simply imposing Af/sol = — ^so\4 = 5 kJ/mol and varying the absolute value of Af/gol.B 

(note that we do not normalize kobs here). An important fact emerges. Whenever the solvation 
free enthalpy is positive the drop of the rate is proportional to exp (| AAGsoi |) (and is hence equal 
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for all curves) because Eq. (??) is basically exact, and the rate is dominated by crossing of the 
polymer gel (Tg Too ^ Tnp). However, when the solvation free enthalpy becomes negative Tg 
can decrease enough to be comparable or even below Too. In this case, the rate is not proportional 
to exp(|AAGsoi|) anymore and since Too has a weaker dependence on temperature than Tg, k^hs be¬ 
comes almost constant. This becomes evident for the curve where Af/soi,s = kJ/mol, where 
the drop in the rate is just around 30% compared to an almost 20-fold decrease for the other cases. 
Finally, panel C) refers to the case where we can control the drop in the diffusion coefficient while 
keeping constant the solvation free enthalpy. In this case, we imposed the same Af/soi^(g) and 
ASsoi,a(s) for all curves (giving AAGgoi = —5 kJ/mol), and we only change the diffusion coeffi¬ 
cient in the B state. A situation similar to that observed in panel A) is seen; The chosen change 
in AGsoi would force an upward jump of the rate at the critical temperature, but this is counterbal¬ 
anced by the drop expected due to a lower diffusion coefficient in the high-temperature, collapsed 
state. Depending on the relative magnitude of these two effects, different signatures are again ob¬ 
served. 

One last thing to notice in Fig.|^is the behavior of the rate as a function of temperature away from 
the transition region. Whenever this is dictated by Tg (i.e. when Eq. (??) is valid), the rate increases 
or decreases with temperature depending on the sign of the solvation enthalpy: a negative enthalpy 
means a decrease with temperature, whereas the opposite is observed for positive values, with the 
rate of decrease (i.e. the slope of the curve) being higher the higher the absolute value of A^fsoi- 
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Figure 5: Top: Normalized reaction rate as a function of temperature, in the case that decreases 
(red curve) or increases (blue curve) upon crossing T^. Note that since the diffusion coefficient in 
the gel is always expected to drop upon polymer collapse at T^, an increase in 1^ can be observed 
only if the solvation free enthalpy decreases, leading to a better interaction of the reacting molecule 
with the solvating gel environment, and hence to an increased concentration of reactants inside the 
gel. AGsoi enters exponentially in the definition of l^(see Eq. (??)), hence even a decrease of a few 
ksT can offset a large reduction in diffusion coefficient. On the contrary, an increase in AGsoi can 
lead to an even larger reduction in the reaction rate than what would be expected by a simple 
reduction in the diffusion coefficient. 


Evidence of the two qualitatively different behaviors we predict based on our theory can be 
found in recent experiments on different PNIPAM-based nanoreactor architectures,^^^^ which we 
show in Fig.|^ In general, one expects that, for a hydrophilic molecule, the solvation free en¬ 
thalpy is lower in the low-temperature, hydrophilic state of the gel compared to that in the high- 
temperature hydrophobic state. Hence, upon crossing Tcfrom below there will be both a drop in 
the diffusion coefficient and an increase in the solvation free enthalpy that will exponentially re¬ 
duce the reaction rate. For a hydrophobic molecule the opposite behavior should be observed for 
the solvation free enthalpy, and an upward jump in the rate can be observed. These cases are in 
full agreement with the experimental data.^i^^ In the first case, Carregal-Romero and cowork¬ 
ers looked at the reduction of hexacyanoferrate(III) (Fe(CN)g“), a strongly hydrophilic molecule 
due to its charge, by borohydride anions using a core-shell nanoreactor.^ These workers indeed 
found a marked drop in the reaction rate (blue squares in Fig.|^, much larger than what would 
be expected due to the drop in the diffusion coefficient as calculated using hydrodynamic theories 
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that should well describe the relevant experimental regime.l^^lS] instead, Wu et o/.used yolk-shell 
nanoreactors to control the reduction of 4-nitrobenzene using borohydride anions (BH^ ). Going 
from below to above the T^, gel permeability is increased in this case due to a lower AGsoi in 
the high-temperature state, inducing the observed jump in the reaction rate (red circles in Fig.|^ 
bottom) despite the decrease in the diffusion coefficient in the collapsed network. An important 
aspect to point out is that the aforementioned reactions are red-ox reactions, and hence in principle 
of second order. However, in the experiments we compare to the concentration of BH 4 was chosen 
in a way to make the reaction pseudo first order (i.e. BH 4 was present in large excess with respect 
to stoichiometric conditions), as also explicitly verified from the experimental kinetic data.^^^^^^ES 
In this case, our theory is still fully applicable. 

Before we conclude, let us now briefly discuss the main limitation of our model. Since we treat 
the reaction within the Debye picture, the polymer cage simply acts as di fixed, external field for the 
diffusing reactant. A more complete description instead would account for the fact that the state 
of the polymer, e.g. its local density, will adapt to the presence of the diffusing specie. Although 
description of the reaction kinetics including such a coupling is possible using techniques such as 
(classical) dynamic density functional theory^ or coarse-grained molecular dynamics simulations, 
the intrinsic complication arising in this more complex case would make the problem analytically 
intractable, and one would need to resort to numerical calculations. Our aim here was instead to 
provide a simple model to rationalize the likely origin of important trends observed in polymer- 
based nanoreactors, as we do. Certainly, a more quantitative description of the problem will require 
the use of such techniques, and we are already planning such simulations for the future. 

CONCLUSIONS 

We developed a model to describe the reaction rate observed in polymer-based, stimuli-responsive 
catalytic nanoreactors. The theory combines a two-state thermodynamic model with the descrip- 
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tion of the reactants’ diffusion which is based on Debye’s theory of diffusion through an energy 
landscape. Model calculations highlight the importance of the solvation free enthalpy difference 
between the bulk solvent and the nanoreactor’s polymeric cage. The theory predicts not only a 
sudden decrease in the observed rate, but also a possibility for rate enhancement, depending on 
the change in solvation free enthalpy at the swollen-to-collapse transition of the PNIPAM-based 
nanoreactors. Such rate enhancement has been observed in recent experiments^i^^ and corrobo¬ 
rates our description. The entire treatment demonstrates that nanoreactors can be used to enhance 
the selectivity of catalysis by nanoparticles. 
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Supplementary Information 

A step-by-step derivation for the full Debye’s equation (including finite reaction rates), specifying 
all assumptions and limitations, is provided in the Supplementary Information. This material is 
available free of charge via the Internet at http://pubs.acs.org. 
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